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ABSTRACT 


The  activities  in  the  reporting  period  focussed  attention 

on  the  following  topics: 

1 . Final  technical  setup  of  the  stationary  lidar  system 
by  improvement  of  the  electronic  data  processing,  full 
automation  of  the  measuring  run  via  a separate  computer 
program  for  tropospheric  and  stratospheric  backscatter 
measurements,  improvement  of  the  photon  counting  system. 

2.  Final  mechanical  and  electronic  setup  of  the  mobile  lidar 
system  including  test  run  of  both  laser  systems,  comple- 
tion through  some  technical  extensions  as  compared  to  the 
original  layout  (e.g.  a heater  for  the  cooling  water 
needed  for  transport  of  the  system  in  winter),  prepara- 
tory works  for  the  inclusion  of  the  primary  neodymium- 
glass  laser  wavelength  (1060  nm). 

3.  Theoretical  and  practical  works  designed  to  derive  lidar 
calibration  constants  for  the  stationary  lidar.  This 
includes  a)  calculation  of  Rayleigh  backscatter  from 
air  density  measured  by  radiosondes  at  very  weak  aerosol 
backscatter  and  b)  in  situ  measurement  of  the  given  size 
distribution  of  the  aerosol  particles  at  our  mountain 
stations  during  simultaneous  lidar  shots,  and  comparison 
with  the  measured  aerosol  backscatter  intensity. 

4.  Numerous  lidar  measuring  series  for  remote  sensing  of 
the  time  variation  of  aerosol  concentration  and  size 
distribution  demonstrating  the  practical  application 
have  been  carried  out  in  the  troposphere  in  the  course 
of  the  day. 

5.  Operation  of  the  fully  developed  single  photon  counting 
system  permitted  the  recording  of  numerous  backscatter 


profiles  (on  the  average  2-3  per  month)  in  the  strato- 
sphere up  to  35  km.  These  measurements  made  it  possible 
to  trace  changes  in  the  aerosol  layer  in  the  lower 
stratosphere  with  consideration  of  time,  height  range, 
cuid  intensity,  and  they  seem  to  indicate  a second 
aerosol  layer  in  the  middle  stratosphere  which  remains 
to  be  verified. 

The  method  used  in  measuring  the  stratospheric  aerosol 
backscatter  can  be  considered  as  satisfactory  and  fully 
developed. 

Work  remains  to  be  done  on  the  calibration  of  the  lidar 
system  since  the  hitherto  existing  data  are  not  yet 
substantiated.  Further  technical  and  theoretical  improve- 
ments turned  out  to  be  necessary.  These  calibration  proce- 
dures are  to  be  finalized  before  realistic  routine  remote 
measurements  of  the  aerosol  parameters  in  the  troposphere 
can  be  started. 


I.  OBJECTIVES  AND  BACKGROUND 


The  background  of  the  investigations  to  be  discussed  in 
the  present  report  can  be  found  in  the  preceding  reports 
R.  REITER  et  al. (1975)  and  R.  REITER  et  al.  (1976  a,  b). 

Our  present  activities  are  a consequent  continuation  of 
the  works  described  in  said  reports. 

The  goals  aspired  to  achieve  in  the  woricing  year  can 
briefly  be  outlined  as  follows: 

1 . Completion  of  the  technical  equipment  of  the  stationary 
(double-frequency  lidar)  system  and  the  on-line  com- 
puter aimed  at  a fully-autoraatic  performance  of  lidar 
shots  and  data  acquisition  concerning  both  tropospheric 
and  stratospheric  (photon  counting)  measurements. 

In  this  connection  elaboration  of  the  necessary  com- 
puter programs. 

2.  Technical  construction  of  the  mobile  lidar  system  (three- 
frequency  lidar)  to  the  most  possible  extent  including 
installation  of  power  supply,  cooling  aggregate,  data 
acquisition  in  the  truck,  and  final  assembly  of  the 
optical  system  (tiltable  transmitter  and  receiver)  on 
the  trailer;  first  measurements. 

3.  Determination  of  calibration  constants  cf  the  lidar  sys- 
tems in  such  a manner  that  through  measurements  of 
aerosol  backscatter  intensity  in  inoi'e  than  one  wave- 
length the  aerosol  particle  concentration  can  be  deter- 
mined simultaneously  in  different  ordereof  magnitude. 

For  the  calibration  we  use  a)  at  very  weak  aerosol  con- 
centration the  Rayleigh  backscatter,  the  intensity  of 
which  can  be  calculated  from  atmospheric  density  di- 
rectly measured  by  means  of  radiosondes,  and  h)  an  in- 
situ  calibration  is  made  of  the  aerosol  backscatter  by 


measurement  of  the  given  aerosol  size  distribution  and 
particle  density  at  our  measuring  stations  at  different 
levels . 

4.  The  aerosol  backscatter  intensity  is  also  to  be  measured 
in  the  stratosphere  by  means  of  photon  counting.  As  the 
same  is  very  weak,  the  Rayleigh  (molecular)  backscatter 
must  be  exactly  known.  For  this  reason  own  rawinsondes 
are  launched  from  which  the  actual  atmospheric  density 
data  can  be  derived. 

5.  Finally,  measurements  as  regularly  as  possible  are  to 
be  performed  also  in  the  troposphere  for  studying  the 
vertical  and  time  variations  of  the  aerosol  concentra- 
tion. For  this  purpose  the  calibrated  lidar  systems 
are  to  be  employed.  As  far  as  the  calibration  proce- 
dures are  not  finished,  such  measurements  will  be  re- 
placed by  the  derivation  and  evaluation  of  relative 
backscatter  profiles. 
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II.  STATUS  OF  THE  LIDAR  TECHNIQUE  AT  THE  INSTITUTE  FOR 
ATMOSPHERIC  ENVIRONMENTAL  RESEARCH 

1 . The  stationary  system 

1.1.  Parts  of  the  system  which  have  not  been  changed  since 

1976 

In  the  last  annual  report  from  July  1976  a complete  and 
detailed  description  of  the  stationary  lidar  system  has 
been  given.  The  following  parts  of  the  system  have  not 
essentially  been  changed  in  the  meantime: 

- the  Q-switched  ruby  laser  transmitter  with  694  nm 
wavelength,  100  MW  peak  output  power  and  1 Hz  max- 
imum pulse  repitition  frequency 

- the  power  supply  for  the  laser  and  the  Q- switch 

- the  fluorescence  suppression 

- the  frequency  doubler  for  generating  the  second  trans- 
mitter wavelength  (347  nm) 

- the  two  energy  monitors  with  peak  digitizer 

- the  fixed  vertical  Cassegrain  receiving  telescope 
with  52  cm  diameter 

- the  EMI  9618  photomultiplier  tube  for  the  detection 
of  the  backscattered  light  (for  medium-range,  analog 
recordings  only) 

- the  Biomation  8100  transient  recorder 

- the  gain  switching  amplifier 

- the  photon  counter  with  range  gating 

- the  preamplifier 

- the  on-line  computer  with  graphic  display  terminal 
and  punch  tape  output. 


1 


It 


These  parts  are  described  in  the  1976  annual  report,  pp,  7 - 


IB,  and  in  the  diagrams  Fig.  III.  1 and  2.  They  are  once 
more  shown  in  tlie  photographs  Figs.  II.  5 - II.  9 of  this 
report. 


1^2.  Improvements  and  extensions  of  the  stationary  system 
aince_July 

During  the  last  year  a number  of  improvements  and  extensions 
to  the  lidar  system  hardware  have  been  performed  which  par- 
tially have  been  announced  in  the  previous  report  (parts  VI 
and  VII).  The  major  part  of  the  innovations  concerns  the 
improvement  of  tlie  data  acquisition  and  the  automation  of 
the  system  operation,  a smaller  part  the  transmitter-receiver 
system. 

Tlie  block  sclieme  Fig.  11.1  presents  the  complete  stationary 
lidar  system  in  its  present  state.  For  the  sake  of  clearness 
not  all  control  lines  are  Included  in  the  figure;  they  are 
shown  completely  in  Fig. II. 2 whicii  will  be  described  below. 
The  innovations  will  now  be  described  in  detail: 

The  FMI  pliotomviltiplj  er  tube  is  not  specially  designed  for 
photon  counting  purposes.  The  quite  low  secondary  emission 
factor  at  the  dyuodes  (2  - 1 , depending  on  the  applied  volt- 
age) causes  heavy  statistical  fluctuations  of  single  photon 
peaks.  The  resulting  broad  pulse  height  distribution  cui've 
renders  an  exact  discrimination  between  signal  and  noise 
pulses  very  difficult  or  impossible  at  all.  So  we  found  a 
slowly  decreasing  background  superposed  on  the  return  sig- 
nals, probably  due  to  the  strong  close-range  signal  in  spite 
of  the  range-gating  of  the  tube,  if  the  full  sensitivity  of 
the  receiver,  without  any  at.tenuatl oia  by  the  gray  wedge,  was 
used.  This  background  covild  only  be  eliminated  by  reducing 
the  received  backsca ttered  light  by  a factor  of  10.  Further- 


more,  the  pulse  frequency  resolution  of  the  EMI  9816  tube 
proved  to  bo  quite  low,  probably  due  to  tlie  high  number  of 
dynodes.  For  these  reasons,  the  acquisition  of  a complete 
backscatter  profile  from  6 to  30  km  altitude  required  the 
application  of  a great  number  of  different  gray  wedge  set- 
tings in  order  to  get  sufficient  counting  statistics  on  the 
one  hand  and  avoid  exessive  pulse  overlapping  on  the  other. 

For  these  reasons  a RCA  8832  photomultiplier  tube  has  been 
provided.  This  tube  is  equipped  with  a gallium  phosphide 
coated  first  dynode,  providing  a secondary  emission  ratio 
of  up  to  50  at  this  dynode  and  consequently  a greatly  im- 
proved counting  statistics.  First  measurements  revealed  the 
performance  of  this  new  tube  to  be  surprisingly  much  better 
than  the  EMI  tube.  Togetlier  with  a new,  faster  jUioton  counter 
input  comparator,  the  pulse  frequency  range  was  extended  by 
a factor  of  no  less  than  10.  Fur  Lher'more , even  with  full  re- 
ceiver sensitivity  no  superposed  background  at  all  is  ob- 
served with  the  return  signals.  Consequently,  the  total  num- 
ber of  laser  shots  required  for  tlie  measurement  of  a complete 
backscatter  profile  with  given  accuracy  is  now  considerably 
reduced . 

Unfortunately,  at  least  at  the  time  of  commission  of  tlie  tube 
no  version  was  available  with  sufficient  sensitivity  ns  well 
for  69^  tim  as  for  3^7  nm  wavoJengtli.  So  for  the  analog  meas- 
urements conducted  in  both  frequencies  the  EMI  tube  is  still 
to  he  used  and  is  to  be  replaced  by  the  RCA  tube  for  strato- 
spheric measurements  performed  with  69/1  nm  wavelengtii  only. 

In  order  to  facilitate  the  change  of  tlie  tubes,  the  h.eavy 
and  unhandy  photomultiplier  cooling  box,  which  proved  to  be 
unnecessary,  was  replaced  by  two  separate  smaller  convention- 
al housings  for  each  tube  with  separate  sockets  and  voltage 
divider  networks  allowing  the  change  of  the  tul’es  within 
less  than  a minute.  Examples  of  stratospheric  measurements 
with  the  new  tube  are  presented  in  part  IV  of  this  report. 
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Pol§rizer_for_meacurement_gf_signal_degolarization 

With  the  initial  single-frequency  version  of  the  system  a 
polarizer  was  available  in  front  of  the  photodetector.  The 
plane  of  polarization  was  changed  by  turning  the  whole  de- 
tector housing.  With  the  onset  of  the  two-frequency  opera- 
tion the  polarizer  was  removed  temporarily  because  the  po- 
larization plane  of  the  3^7  nm  laser  beam  is  perpendicular 
to  that  of  the  69^  nm  beam,  and  the  change  of  the  transmit- 
ter frequency  would  have  required  each  time  the  time-con- 
suming turning  of  the  photomultiplier  case  by  90  degrees. 
However,  as  the  measurement  of  the  signal  depolarization  is 
desirable  for  several  reasons,  a motor-operated  turnable 
mount  for  the  polarizer  was  fixed  to  the  receiver,  which 
may  be  controlled  from  inside  either  manual  or  by  the  comput- 
er (see  below). 

1 |§9lSgCgyDd_stog_tube_fgr_the_receiver 

The  receiver  sensitivity  during  daytime  is  limited,  among 
other  things,  by  the  skylight  background  level.  Initially 
skylight  bypassing  the  telescope  secondary  mirror  and  en- 
tering directly  the  photomultiplier  cathode  contributed 
very  essentially  to  the  background  level.  For  this  reason 
a tube  T was  mounted  at  the  center  hole  of  the  primary  mir- 
ror (see  Fig.  II.  1)  which  reduced  the  background  by  a fac- 
tor of  nearly  five.  This  resulted  in  a marked  improvement 
of  accuracy  especially  of  the  weaker  347  nm  signals. 

li2i4i__Autgmatic_ temperature _contrgl_for_the_laser 

The  laser  is  cooled  by  a closed  water  circuit  which  is  re- 
cooled by  tap  water  by  means  of  a heat  exchanger.  The  exact 
wavelength  of  the  ruby  laser  depends  somewhat  on  temperature. 
In  order  to  avoid  unintentional  tuning  of  the  laser  to  one 
of  the  near-by  water  vapor  absorption  lines,  a magnetic  valve 
operated  by  a contact  thermometer  was  applied  to  the  tap  wa- 
ter line.  This  valve  stops  the  tap  water  flow  if  the  temper- 
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ature  of  the  secondary  circuit  water  drops  below  a preset 
value. 

EloEBY_^i§9_5§>DorY 

In  addition  to  the  16  K core  memory  a floppy  disc  memory  is 
now  available  with  the  on-line  computer.  It  consists  of  two 
separate  units  with  300  k capacity  each  and  avoids  the 
necessity  of  intermediate  punch  tape  output  of  averaged 
signals  during  a series  of  measurements  (e.g.  between  a 
694  nm  and  a 347  nm  measurement).  Instead  of  that,  any  num- 
ber of  single  signals  may  be  stored  on  one  or  more  floppy 
discs  and  output  on  punch  tape  later  after  finishing  a com- 
plete series  of  measurements.  This  device  is  also  very  con- 
venient for  storage  of  programs.  The  usage  is  described 
more  detailed  in  part  II. 1.3.  of  this  report. 

liii§i__^|y!tral_filter_wheel 

Previously  the  receiver  sensitivity  attenuation  was  accom- 
plished by  a continuous  neutral  or  gray  wedge,  remotely 
operated  by  a magslip.  Because  it  proved  to  be  difficult 
setting  the  wedge  to  a certain  attenuation  value  with  suf- 
ficient reproducibility,  it  has  been  replaced  by  a motor- 
operated  filter  wheel  with  3 discrete  neutral  filters  and 
one  empty  hole  (zero-attenuation  position). 

1 06Yi9®§_fQE_Q9n?Ei®t§_iy^2!B§Jio_o2eration_of  the 

iy|tera 

The  Biomation  8100  transient  recorder  is  fully  externally 
programmable.  For  utilization  of  this  possibility  a special 
computer  interface  has  been  available  for  more  than  a year. 
It  has  already  been  mentioned  in  the  July  1976  annual  re- 
port and  is  shown  in  Fig.  III.1,  denoted  as  "Interface  2". 
However,  the  external  programming  of  the  transient  recorder 
only  seemed  quite  unsatisfactory,  and  so  it  was  decided  to 
develop  additional  electronic  circuits  allowing  a complete 
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automatic,  computer-controlled  operation  of  the  whole  system. 
This  should  include  the  following  single  operations: 

- programming  of  the  transient  recorder 

- laser  start 

- change  of  transmitter  wavelength 

- change  of  interference  filter 

- change  of  neutral  filter 

- setting  of  photomultiplier  high  voltage 

- setting  of  range  gate  length  and  delay 

- triggering  of  the  data  acquisition  electronics 
for  the  measurement  of  the  skylight  background 

- change  of  the  polarizer  position. 

The  control  electronics  which  has  been  developed  here  at  the 
institute  by  our  co-workers  (see  Fig.  II. 2)  will  now  be  des- 
cribed in  detail,  including  a modified  data  transfer  from 
transient  recorder  and  photon  counter,  guided  by  Fig.  II. 2. 

liiiZiIi==i22y|Zoy|Ey|_i2S|£2i_||cilitie|_of_the_comput|r 

For  the  exchange  of  data  between  computer  and  external  devices 
two  32-bit  input/output  controllers  are  available.  The  input 
and  output  lines  are  divided  into  four  groups  or  channels 
d 8 bit  each,  which  are  accessed  by  channel  select  function 
codes  from  the  computer.  During  an  output  transfer,  the  8-bit 
word  from  the  computer  is  loaded  into  the  input  register  of 
the  selected  channel,  and  an  output  data  strobe  pulse  is 
issued  to  the  device.  During  an  input  transfer,  the  8 input 
lines  from  the  selected  channel  are  scanned  by  an  input 
multiplexer  and  the  resulting  input  word  is  supplied  to  the 
computer  when  "data  in"  occurs.  Additonally  an  external  in- 
terrupt line  is  available  with  each  I/O  controller. 

li2iZi|i__|^|0iii2^=£f22ri®r_2£gg£i5iiDS_§24_data_transfer 

All  front  panel  controls  of  the  Biomation  8100  transient  re- 
corder may  be  remotely  controlled  without  regard  to  the  phys- 
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ical  position  of  the  manual  buttons  and  knobs.  The  digital 
control  is  effected  via  bit-parallel,  word-serial  data  ex- 
change. The  control  functions  of  the  transient  recorder  are 
divided  into  several  groups,  such  as  arm  or  trigger  control 
groups,  to  which  the  remote  control  can  be  applied  separate- 
ly. Thus  any  control  group  or  combination  of  control  groups 
may  be  under  remote  digital  control,  other  groups  may  be  re- 
leased for  manual  front  panel  operation. 


The  remote  control  of  the  recorder  requires  1 6-bit  control 
words,  which  are  supplied  by  interface  No.  2 from  the 
computer.  By  means  of  this  interface  the  1 6-bit  word  is 
generated  from  two  8-bit  words  issued  sequentially  from  the 
computer  to  the  interface.  This  is  shown  in  Fig.  II. 2,  where 
heavy  lines  indicate  a set  of  more  than  one,  usually  8 or  9, 
individual  data  lines. 


I 

\ 
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The  digital  data  stored  in  the  transient  recorder  are  out- 
put sequentially  as  8-bit  words.  Previously,  the  data  of 
both  transient  recorder  and  photon  counter/peak  digitizer 
had  been  transferred  to  the  computer  via  interface  No.  1, 
as  described  in  detail  in  the  July  '76  annual  report.  With 
the  computer  a single  channel,  8-bit  interface  was  used  for 
this  purpose  which  otherwise  is  destined  for  punch  tape  in- 
put and  output  . As  with  the  I/O  controllers  sufficient  input 
channels  are  available,  one  channel  of  controller  No.  1 is 
now  used  for  the  data  transfer  from  transient  recorder  and 
photon  counter.  Data  from  photon  counter/photodiode  peak 
digitizer  are  still  transferred  via  interface  1 , the  tran- 
sient recorder  data  via  the  new  interface  2.  The  usage  of 
this  interface  offers  several  advantages.  First,  data  re- 
cording may  be  done  in  the  "off"  mode  of  the  recorder,  al- 
lowing the  examination  of  the  data  prior  to  the  transfer  to 
the  computer.  The  data  transfer  is  initiated  by  the  computer 
by  issuing  a special  instruction  word  to  the  recorder. 
Furthermore,  for  technical  reasons  the  data  transfer  is 
much  faster  than  with  the  former  interface,  which  allows 
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the  transfer  routine  for  one  data  word  to  be  executed  by 
the  computer  within  less  than  100  microseconds.  Consequent- 
ly, as  the  recorder  memory  made  up  of  MOS  shift  registers 
needs  refreshing  between  output  words  at  read  rates  slower 
than  100  microseconds  per  word,  the  time  required  for  the 
transfer  of  the  total  memory  contents  is  reduced  from  2.5 
to  about  .2  seconds. 

The  usage  of  only  one  input  channel  for  transferring  both 
the  transient  recorder  and  photon  counter  data  is  made 
possible  by  tri-state  outputs  in  both  interfaces.  An  inter- 
nal electronic  switch  operated  by  a computer  instruction, 
which  is  not  included  in  Fig.  11.2,  prevents  simultaneous 
output  from  both  interfaces. 

Rr^viously  the  range  gate  length  ("gate”)  and  delay  time 
("delay")  had  been  set  manually  by  two  4-digit  decimal 
switches  each,  in  increments  of  .1684  microseconds.  The  8 
digits  may  now  be  set  by  the  computer,  where  for  saving  of 
output  channels  multiplex  operation  is  used.  The  most  sig- 
nificant 4 bits  of  the  8-bit  word  output  from  the  computer 
are  used  as  digit  address,  the  4 least  significant  bits 
denote  the  "value"  of  the  corresponding  digit,  BCD-coded. 
Manual  control  of  the  gate  and  delay  is  still  possible. 

The  change  from  manual  to  remote  operation  is  accomplished 
either  via  the  computer  again  or  via  an  external  switch, 
which  is  not  included  in  the  block  scheme. 

1^2i2i^f  i [iigb_y2iiage_control 

The  output  voltage  of  the  Ortec  photomultiplier  high  voltage 
supply  may  be  controlled  externally  by  applying  a DC  refer- 
ence voltage  between  0 and  9 Volts  to  a special  connector. 
This  voltage  is  generated  by  an  analog  multiplexer,  which 
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acts  like  a manual  step  switch  on  a precise  voltage  divider 
chain.  The  analog  multiplexer  is  controlled  by  the  computer 
via  the  above-mentioned  multiplexer.  Only  hundreds  and  thou- 
sands of  volts  can  be  set  remotely,  which  has  been  decided 
to  be  sufficient.  Manual  remote  control  of  the  high  voltage 
is  possible,  similarly  as  with  the  range  gate  control.  Two 
digits  of  a four-digit  decimal  switch  are  used  for  this  pur- 
pose (cf.  the  photograph  Fig.  II.9).  Precautions  are  made  to 
avoid  the  unintentional  application  of  a too  high  voltage  to 
the  photomultiplier  tube.  If  an  attempt  is  made  to  set  the 
voltage  to  a value  above  a limit  voltage  either  manually  or 
by  the  computer,  the  voltage  is  set  to  zero  and  an  error 
message  is  issued  to  the  computer  via  a combined  error  line. 
Error  signals  from  other  devices  are  also  transferred  through 
this  line.  The  limit  voltage  may  be  changed  by  a computer- 
operated switch  from  1300  Volts  for  analog  measuring  to 
2000  Volts  for  photon  counting.  After  any  change  of  the  HV 
a certain  settling  time,  depending  on  the  amount  of  the 
change,  is  required  by  the  PMT  to  attain  its  nominal  gain. 

In  order  to  prevent  measurements  before  the  gain  has  stabi- 
lized an  error  signal  is  output  for  a time  interval  simu- 
lated by  a R-C  network. 


liiiZi5i__N6y^C§i_fii£§P_wheel_control 

As  mentioned  above,  the  previously  used  continuous  neutral 
wedge  has  been  replaced  by  a filter  wheel  containing  3 neu- 
tral filters  with  7.7,  1.5,  and  .3  % transmissivity,  respec- 
tively, and  one  empty  hole  ("no  attenuation"  position  of  the 
wheel).  The  wheel  is  operated  by  a servo  motor.  The  position 
of  its  axis  is  indicated  by  a potentiometer  connected  with 
it.  By  applying  a logical  "1"  signal  to  one  of  the  four  in- 
put lines  (denoted  by  0...3  in  Fig.  II. 2)  of  the  control 
unit  by  the  computer,  the  servo  motor  moves  the  wheel  to  the 
desired  position.  An  error  signal  is  output  as  long  as  the 
motor  is  moving.  Manual  operation  is  possible  by  means  of  a 
toggle  switch  with  positions  0-3  and  "ext".  In  the  latter 
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position  only  remote  control  is  possible. 
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1 Xransmitter_freguency_and_golarizer_cgntrol 

As  the  polarization  planes  of  the  694  and  347  nm  laser  beams 
are  perpendicular  to  each  other,  the  0°-position  of  the  po- 
larizer for  694  nm  is  physically  the  same  as  the  90°-position 
for  347  nm,  with  the  angle  denoting  that  between  the  planes 
of  the  polarizer  and  of  the  polarization  of  the  incident 
light.  As  it  was  desired  to  maintain  the  position  of  the  po- 
larizer with  respect  to  the  polarization  plane  of  the  trans- 
mitted light  when  the  transmitter  wavelength  is  changed,  a 
combined  control  unit  has  been  developed.  A wavelength 
change  thus  requires  three  single  operations;  Shifting  of 
the  frequency  doubler  into  the  laser  beam  path,  changing 
the  interference  filter  and  changing  the  polarizer  position. 
The  frequency  doubler  servo  motor  has  hitherto  been  operated 
by  a two-position  switch  at  the  operating  console;  the  end 
positions  of  the  doubler  carriage  have  been  indicated  by 
two  lamps  operated  by  micro  switches.  The  control  has  now 
been  removed  from  the  console  switch  to  a computer-operated 
relay  in  the  control  unit.  When  the  carriage  is  moving,  i.e. 
when  none  of  the  micro  switches  is  closed,  an  error  signal 
is  output  to  the  computer. 

The  interference  filter  wheel  has  formerly  been  moved  by 
means  of  magslips.  As  the  magnetic  momentiom  goes  to  zero 
when  the  wheel  approaches  the  desired  position,  the  posi- 
tioning is  not  precisely  reproducible.  (As  mentioned  above, 
the  magslip  drive  for  the  neutral  wedge  had  been  replaced 
by  a servo  motor  for  the  same  reason).  The  magslip  has  been 
replaced,  therefore,  by  a rotary  magnet.  This  device  is  oper- 
ated from  the  control  unit. 

The  polarizer  drive  is  quite  similar  to  that  of  the  frequency 
doubler  carriage.  It  is  moved  by  a DC  servo  motor,  which  is 
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stopped  by  microswitches  as  soon  as  one  of  the  two  end  posi-  <i 

tions  is  reached.  The  motor  is  operated  by  the  frequency/  \ 

polarizer  control  unit  which  is  in  turn  controlled  from  the  j 

computer  via  four  lipes  which  are  denoted  in  Fig.  II. 2 by  j 

"R%  "UV",  "O®",  and  "90°”.  Applying  a logic  "one"  signal  I 

to  the  "R"  line,  which  means  694  nm  (red),  the  frequency 
doubler  is  moved  out  of  the  laser  beam  path,  the  694  nm  in- 
terference filter  is  moved  into  the  telescope  focus  and  the  ‘f 

polarizer  either  parallel  to  the  polarization  plane  of  the  ■; 

694  nm  incident  signal  if  a "one"  is  applied  to  the  "0°"  J 

line,  or  normal  to  it  if  a "one"  is  issued  to  the  "90°"  line.  ;| 

Accordingly,  if  a "one"  is  applied  to  the  "UV"  (ultraviolet)  | 

line,  the  frequency  doubler  is  moved  into  the  laser  beam,  the  | 

347  nm  filter  into  the  focus,  and  the  polarizer  parallel  to  j 

the  347  nm  plane  if  a "one"  is  issued  to  the  "0°"  line,  or  || 

normal  to  the  347  plane,  if  a "one"  is  issued  to  the  "90°"  i 

line.  The  actual  position  of  the  polarizer  depends,  there-  \ 

fore,  on  both  the  polarizer  and  the  frequency  control  signal  j 

from  the  computer.  ] 

J , 

The  four  control  signals  may  also  be  issued  to  the  control  \ 

unit  manually  via  pushbuttons.  \ 

The  electronics  for  remote  control  of  range  gate,  high  ] 

voltage,  filter  wheels,  polarizer  and  frequency  doubler  is  i! 

contained  in  a 19"  chassis.  Its  front  panel  is  shown  in  the 

/I 

photograph  Fig.  I I. 9.  i 


{ 

/' 
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The  laser  is  started  by  closing  a normal-open  contact  at  the 
laser  power  supply  (cf.  description  in  the  previous  reports).  [■ 

This  contact  may  now  be  closed  remotely  by  a relay  operated  'i 

from  the  computer  via  interface  No.  2.  / 

f 

i' 

1.2. 7. 8.  External  clock  J 

= = = sss  = = = =:  = =;  = = s = sss  = sss=;  ^ 

An  external  digital  clock,  synchronized  by  line  frequency,  j 
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has  been  developed.  The  four  digits  (hours  and  minutes)  are 
output  in  BCD  code  and  input  to  the  computer.  Two  8-bit 
channels  of  one  I/O  controller  are  used  for  this  purpose. 
The  clock  may  be  questioned  by  the  computer  and  the  actual 
time  output  on  punch  tape  together  with  the  signals,  as 
described  in  detail  in  paragraph  1.3.  Furthermore,  measure- 
ments may  be  started  automatically  at  any  time  by  means  of 
the  clock. 


An  easily  accessible  "emergency"  pushbutton  has  been  in- 
stalled which  gives  an  interrupt  signal  to  the  external 
interrupt  input  of  one  I/O  controller.  Pressing  this  but- 
ton interrupts  the  current  program  if  the  computer  inter- 
rupt system  is  enabled  (see  below).  This  is  used  for  ex- 
ample when  a laser  shot  is  to  be  repeated  because  of 
superposed  spurious  pulses. 


liii.Zil2i__^y^2iy5ic_trigger_and_reset_of_phgtqn_cgunter 

The  automatic  background  measurement  with  the  photon  counter 
as  well  as  with  the  transient  recorder  requires  the  computer 
triggering  of  the  photon  counter  (in  the  latter  case  because 
of  the  needed  range  gate  start).  This  is  now  possible  via 
one  of  the  data  output  lines  from  the  computer.  Simultane- 
ously with  the  photon  counter  triggering  a pulse  is  fed 
through  a cable  to  one  of  the  external  trigger  input  con- 
nectors of  the  transient  recorder. 

Via  another  data  output  line  the  combined  photon  counter/ 
peak  digitizer  unit  may  be  reset  externally.  This  is  con- 
venient for  analog  measurements  as  in  this  case  the  energy 
monitor  data  only  are  needed  and  the  data  transfer  may  be 
stopped  by  reset  when  the  peak  digitizer  readout  is  com- 
pleted. 


> 
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1 . Ti}®_D§^_£2?Byier_programs 

With  the  availability  of  the  automatic  control  of  the  system 
the  development  of  completely  new  programs  was  necessary. 

As  now  by  means  of  the  floppy  disc  memory  programs  can  be 
loaded  very  quickly  into  the  computer,  two  separate  programs 
for  analog  and  photon  counting  measurements  were  prepared. 
They  will  now  be  described  briefly  with  the  aid  of  Figs.  II. 
3 and  4. 


The  flow  diagram  Fig.  II. 3 a and  b shows  the  program  for  au- 
tomatic analog  measurements  with  the  transient  recorder.  The 
automatic  operation  of  the  system  requires  a set  of  control 
words  for  the  transient  recorder  and  the  control  devices. 
They  are  stored  in  a table  which  is  included  in  the  program 
and  loaded  into  computer  together  with  it.  As  like  previous- 
ly one  complete  measurement  series  actually  consists  of 
four  single  measurements,  e.g.  one  low  resolution  measure- 
ment (with  low  photomultiplier  gain)  for  the  lower  parts 
of  the  total  altitude  interval  of  interest  and  one  high- 
resolution  measurement  for  the  upper  parts  in  both  trans- 
mitter wavelengths,  the  table  of  control  words  is  divided 
into  four  groups  of  partially  different  words  for  the  same 
control  operations.  The  table  contains  a standard  set  of 
control  words  suitable  for  most  atmospheric  conditions,  but 
it  may  deliberately  be  changed  from  the  keyboard,  for  ex- 
ample if  extremely  intense  backscatter  requires  a reduced 
sensitivity  of  the  recorder,  if  the  altitude  range  is  to  be 
changed  and  so  on.  Furthermore  it  is  possible  to  store  any 
number  of  different  tables  on  a floppy  disc  and  to  load  one 
of  them  into  the  computer  instead  of  the  original  table. 

Before  starting  the  program  the  actual  date  is  entered  into 
the  computer  from  the  keyboard  for  later  output  on  punch 
tape  and  display  on  the  terminal  screen.  The  program  starts 


with  a decision  whether  a manual  start  of  the  measurements 
is  desired  or  if  an  automatic  start  at  certain  times,  given 
by  the  external  clock,  is  desired.  This  decision  is  done  ac- 
cording to  the  position  of  one  of  the  four  sense  switches  of 
the  computer.  In  the  case  of  automatic  start,  the  computer 
checks  the  current  time  and  compares  it  with  the  desired 
starting  times,  e.g.  every  half  and  full  hour,  and  goes  on 
if  this  condition  is  met.  The  program  then  takes  the  first 
from  the  four  sets  of  control  words  and  stores  them  into  a 
"current  control  word  buffer".  The  subroutine  called  "meas- 
ure", which  the  computer  executes  subsequently,  takes  the 
control  word  from  this  buffer. 

After  execution  of  this  subroutine  which  is  described  below, 
the  program  continues  storing  the  second  set  of  control  words 
into  that  buffer,  again  executes  the  subroutine  and  so  on 
until  all  four  single  measurements  are  completed.  Finally 
the  transient  recorder  and  control  devices  are  returned  to 
front  panel  (manual)  operation  and  the  program  goes  to  start 
location  again. 

The  subroutine,  which  really  represents  the  major  part  of 
the  complete  program,  is  shown  in  Fig.  II. 3 b.  At  first  the 
control  words  from  the  current  buffer  are  output  to  the  de- 
vices, and  all  buffers  for  data  to  be  entered  as  well  as  the 
buffer  containing  the  current  number  of  laser  shots  are  set 
to  zero.  As  soon  as  a possible  error  signal  is  reset,  i.e. 
if  all  servo  motors  are  at  their  correct  positions  and  the 
PMT  high  voltage  has  settled,  the  data  acquisition  elec- 
tronics is  triggered  for  background  measurement.  The  stored 
data  is  then  transferred  to  the  computer.  The  data  transfer 
routines  for  transient  recorder  and  photon  counter/peak  dig- 
itizer are  quite  similar  to  those  in  the  former  program  des- 
cribed in  the  1976  annual  report  and  are  therefore  not  in- 
cluded in  the  diagram.  From  a part  of  the  transient  recorder 
data  the  background  level  is  calculated  and  added  to  the 
background  buffer.  The  next  step  is  firing  the  laser.  There- 
after, the  stored  signal  may  be  viewed  on  the  screen  of  the 


recorder  display  during  an  1 sec.  pause.  The  interrupt  sys- 
tem is  enabled  during  this  pause,  and  by  pressing  the 
"emergency"  button  the  last  laser  shot  may  be  repeated 
(also  the  foregoing  background  measurement)  if  necessary. 

When  the  pause  is  over  the  stored  data  are  transferred  from 
the  transient  recorder  to  the  computer  and  added  to  the  sig- 
nal buffer  contents.  Likewise  the  energy  monitor  value  is 
transferred  and  added  to  the  peak  digitiser  buffer  consents. 
This  is  followed  by  incrementing  the  shot  number  and  com- 
paring it  with  the  intended  number  of  shots.  If  the  current 
number  of  shots  is  still  less  than  the  proposed  number,  the 
background  and  signal  measurements  are  repeated;  otherwise, 
the  time  is  again  input  into  the  computer,  the  background 

sum  subtracted  from  the  signal  sum,  the  "net"  signal  multi- 
2 

plied  by  R , divided  by  the  average  laser  enei’gy,  multiplied 
by  a scale  factor  and  plotted  on  the  terminal  screen  together 
with  the  correct  altitude  scale.  The  current  date  and  time, 
transmitter  wavelength,  range  gate  data  and  transient  recor- 
der input  range  and  trigger  delay  are  also  displayed  and  the 
complete  plot  hardcopied  if  wanted.  The  routine  is  then 
finished  by  storage  of  the  net  signal  on  a floppy  disc  for 
later  punch  tape  output. 


1 

The  flow  diagram  Fig.  II. 4 presents  the  new  photon  counting 
program.  Initially  the  photomultiplier  HV  and  its  limit  is 
set  to  2000  Volts.  Then  the  first  set  of  control  v/ords  for 
range  gate  neutral  filter  wheel  and  polarizer  position  is 
output  from  the  computer.  All  control  words  are  contained 
in  a table  similar  to  the  analog  measuring  program.  After 
setting  the  current  shot  number  and  other  buffers  to  zero, 
the  measurement  is  started  by  triggering  the  photon  counter 
without  firing  the  laser  for  background  measurement.  Since 
the  computer  triggering  is  possible  we  measure  the  back- 
ground in  this  way  instead  of  using  the  eleventh  or  i-efer- 
ence  channel  of  the  counter  which  results  in  an  increased 
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mnl  nuinlu'r;!  and  it  i ;!i!]  .ayi'd  on  tin'  tormiinil  ;’.croon.  Tito  not 
coutit;!  ant'  .aditi'd  tii  ;!  I pn.a  1 buft'or  contont,;'.  (which  i ;;  coi'o 
l;i1ti;illy)  and,  with  I ho  t'xi't'plion  of  tlu'  fin;!t  ;!hot  , t lu' 

;i  t a t I ;i  t 1 ca  1 oi’i'oi';!  of  t ho  .accumu  1 !i  t t'd  count;-!  .ai'o  I cu  1 ;i  t t'd  . 
'I'hon  tho  iK'xt  in«'a:!ui'oinon  t cycle  i ;!  ttt.antt'd  with  tin'  h.ack- 
pnound  mo;!;!Ui'i'mont  . Tho  ;-!!'ni('.'!  of  moat’.unonn'nt ;!  with  a pivt'n 
nanpo  pato  and  lu'uti'.al  filter  ;-!('ttinp  ia  ;!toppod  automatical- 
ly if  the  .avi'i’ap!'  ;!  t .at  i ;■!  t 1 c.a  1 ('I'i’or  i :!  below  h %,  if  *'0  .■!inplo 
;!hot;i  h.avo  hc'on  fired,  or  m.anua  1 1 y by  mi'.an;’.  of  tin'  fourth 
I'ompuloi’  ;!on;n'  .twitt-h.  Tin'  ;ihot  numht'r,  ;!um  of  onci’py  mom - 
lor  value:!,  I'anpo  p.ilt'  d;!t.a  .and  IniU'  .ai't'  tln'ii  outinit  on  I'uin'h 
tape,  .'ind  tin'  rantn'-coria'c  t t'd  .and  ('tn'rr.y-uoi'ma  I i -.n'd  \'.!!'ti;vl 
h.acU;n'.a  I I ('!•  pia'I'ili'  dl:!pl.av('d  i-'ii  I ho  tcianinal  in  .a  :n'mi-lcip 
;!C.'ii(',  lopi'lhcr  with  I in'  I'l’niu'  limit;!. 

Tin'  comploli'  ;;i'i'ii';!  of  mo;i!!ui-t'iiU'n  I 1;!  now  la'pt'alt'd  with  tin' 
folK-iwi’ip  ;!ol;!  I'f  i-onlriil  will'd;!  !inlil  tho  end  of  t hi'  l;itilo 
i ;!  ro.achod.  A;!  in  the  c.a;!c  of  .an.alop  mo.a.'nii' i np,  tlio  fin;il 
evaluation  of  the  d.at.a  i ;!  :!l  i 1 1 porl’oi'mod  off-line  on  the 
mail!  computi'i'  with  \'roci:!o  ploltinp  fac  i I i I i oi! . 

In  peneral  wo  uae  a hoipht  ro:!Olulion  of  ('00  m with  tho 
photon  count  inp  me.aauromon  t,;! , The  corro;;ponii  I np  A mict'o- 
aocond;!  pale  lenplh  per  countinp  channel  i a included,  there- 
fore, in  the  ;!l.anii.ird  control  word  table,  hut  m.ay  In'  chanped 
do  1 i hi'i'.a  I o 1 V via  the  koyho.ard  terminal.  The  at  and.ard  alli- 
I udo  interval;!  ;inil  neutral  filter  aetlinp;!  for  tho  p.arti.al 
meaaurlnp  ;ioi'ioa  are  tho  followinp: 
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Series  Altitude  Interval  Neutral  Filter 

No.  (km)  Attenuation 


1 

34.8 

- 

40.8 

0 

% 

2 

30.0 

- 

36.0 

0 

% 

3 

25.2 

- 

31.2 

0 

% 

A 

20.4 

- 

26.4 

0 

% 

5 

15.6 

- 

21 .6 

92.3 

% 

6 

10.8 

- 

16.8 

92.3 

% 

7 

6 . 0 

- 

12.0 

98.5 

% 

8 

32.4 

- 

38.4 

0 

% 

9 

27.6 

- 

33.6 

0 

% 

10 

22.8 

- 

28.8 

0 

% 

1 1 

18.0 

- 

24.0 

0 

% 

12 

13.2 

- 

19.2 

92.3 

% 

13 

8.4 

- 

14.4 

92.3 

% 

14 

3.6 

- 

9.6 

98.5 

% 

As  shown  in  the  table,  the  complete  series  of  amplitude  inter- 
vals consists  of  two  parts  which  are  shifted  by  300  m with 
respect  to  each  other  resulting  in  a real  height  resolution 
of  300  instead  of  600  m.  Examples  of  backscatter  profiles 
measured  by  means  of  the  photon  counting  technique  will  be 
presented  in  part  IV. 2 of  this  report. 


2.  The  mobile  lidar  system 

The  mobile  lidar  system  has  been  described  in  detail  in  the 
1976  annual  report.  In  the  meantime  the  electronics  and 
auxiliary  equipment  as  cooling  system  etc,  has  been  fixed  to 
the  van.  For  safety  reasons  this  work  had  to  be  dojie  very 
carefully  and  took  several  montlis.  To  prevent  fi'ee?'.ing  of 
the  water  in  the  cooling  circuit  during  transport  of  ttie  sys- 
tem at  air  temperatures  below  zero,  a heater  has  Ix'eti  con- 
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structed  to  which  the  cooling  water  tubings  can  be  connected 
and  which  contains  a little  pump  to  circulate  the  water  slow- 
ly. The  heating  power  is  delivered  by  the  vehicle  battery. 

As  well  as  the  lasers  and  the  electronics  have  been  tested 
separately,  some  defects  especially  in  the  electronics  have 
been  eliminated.  First  testing  measurements  will  be  possible 
as  soon  as  the  final  alignments  are  completed. 

The  mobile  lidar  system  is  shown  on  the  photo  Fig.  11.10. 

It  has  been  presented  at  the  "laser  77"  exhibition  in 
June  at  Munich. 


III.  ABSOLUTE  CALIBRATION  OF  THE  LIDAR  SYSTEM 


1 . Calibration  constants 

The  method  of  absolute  calibration  of  the  monostatic  double 
frequency  ruby  lidar  system  has  been  described  in  previous 
reports  by  REITER  et  al.  (1975,  1976a).  Several  new  cali- 
bration experiments  have  been  carried  out  in  order  to  de- 
termine the  lidar  system  constants  needed  for  the  evalua- 
tion of  lidar  backscattering  profiles  with  respect  to 
aerosol  concentration  and  size  distribution. 

Table  III.1  contains  a list  of  lidar  system  constants 
that  have  been  found  until  now  (see  also  Fig.  III.1). 
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Table  III.1 


Date  A(694  nm)"'^  A(347  nm)^^  fitted  to 


05  Mar 

76 

3.4 

10^^ 

7.7 

10^0 

Wank 

05  Mar 

76 

1 .8 

lo”^ 

6.5 

lo”^ 

Rayleigh 

1 1 Mar 

76 

2.1 

10^^ 

5.1 

10*^ 

Wank 

01  Jul 

76 

2.3 

10^ 

4.2 

io'° 

Wank 

01  Jul 

76 

2.2 

lo'’^ 

3.5 

lo'^ 

Rayleigh 

08  Sep 

76 

2.8 

10^^ 

1.4 

10^“ 

Wank 

04  Mar 

77 

1.4 

10^^ 

2.1 

10™ 

Wank 

09  Mar 

77 

2.0 

lo”^^ 

2.0 

10™ 

Wank 

22  Mar 

77 

2.6 

lo”’'^ 

1.2 

10™ 

Rayleigh 

20  May 

77 

2.9 

10^^ 

8.5 

10™ 

Wank 

1 ) -13 

^ The  unit  of  the  system  constant  is  mV  J m sterad. 

2 ) 

' The  lidar  backscattering  profiles  were  fitted  to  ex- 
pected backscattering  values  computed  from  additional 
data: 

Wank  : combined  Rayleigh  and  aerosol  scattering 
based  on  impactor  measurements  of  the 
aerosol  size  distribution  at  Wank  peak 
(1780  m), 

Rayleigh:  pure  Rayleigh  scattering  deduced  from 
simultaneous  radiosonde  data  during 
days  with  extreme  weak  aerosol  back- 
scatter. 

There  is  a considerable  spreading  in  the  calibration  con- 
stants of  about  40  % at  694  nm  and  more  than  one  order  of 
magnitude  at  347  nm.  Of  course,  the  larger  spreading  in  the 
ultraviolet  is  mainly  due  to  the  higher  fraction  of  Ray- 
leigh scattering.  But  there  are  still  other  uncertainty 
factors  that  are  not  yet  eliminated  by  data  normalization 
and  the  fitting  procedure.  These  uncertainty  factors  may 
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be  of  technical  origin,  or  they  may  enter  the  calibration 
routine  through  the  various  models  and  assumptions  that 
have  to  be  made. 

On  22  March  1977  the  system  constants  were  derived  under 
rather  clean  air  conditions.  According  to  the  impactor 
measurements,  the  aerosol  concentration  was  less  than  10  96 
of  the  "normal"  concentration  measured  at  Garmisch-Parten- 
kirchen  and  at  Wank  peak.  Under  these  conditions  Rayleigh 
scattering  is  the  predominant  scattering  process,  which  is 
described  very  well  by  theory.  Therefore  the  resulting 
lidar  calibration  is  not  so  much  affected  by  assumptions 
on  aerosol  scattering  but  rather  by  apparative  factors. 

The  extremely  high  system  constant  at  3^7  nm  on  this  spe- 
cial day  was  probably  due  to  the  fact  that  energy  moni- 
toring via  photodiode  gave  too  small  transmitted  energies. 
Energy  normalization  then  leads  to  an  additional  factor 
contributing  to  the  lidar  calibration  constant.  It  is 
planned  to  improve  the  present  energy  monitoring  system. 

On  4 March  1977  a relatively  small  system  constant  at  694 
nm  is  found.  The  meteorological  parameters  and  a charac- 
teristic visibility  indicated  the  advection  of  Sahara  dust. 
This  was  confirmed  by  aerosol  samples  on  filters  exposed 
at  the  mountain  stations  and  by  impactor  measurements.  The 
size  distribution  (Fig.  III. 2)  shows  an  enhancement  of  par- 
ticles in  the  size  range  0.4  - 1.0  ^um  diameter.  The  bimodal 
log-normal  function  which  we  use  to  describe  aerosol  size 
distributions  has  a gap  just  in  this  size  range  (Fig.  III. 3). 
Since  backscattering  at  both  694  and  347  nm  is  sensitive 
to  particles  of  that  "medium"  size,  any  deviations  from  the 
measured  particle  size  spectrum  will  alter  the  calibration 
constants  considerably.  Moreover,  the  shape  of  these  par- 
ticles, their  refractive  index  and  its  humidity  dependence 
are  probably  different  from  the  mean  aerosol  to  which  our 
model  refers.  A variation  of  these  parameters  is  not  yet 
possible  with  this  calibration  routine. 
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Fig.  III. 3 shows  a bimodal  log-normal  distribution  with  the 
same  backscattering  properties  at  69^  and  3^7  nm  as  the  meas- 
ured size  distribution  (solid  line)  measured  on  4 Mar  77. 
Since  we  use  log-normal  distributions  with  fixed  mean  dia- 
meters (0.4  and  2.0  ^um)  and  standard  deviation  0.3,  a 
size  distribution  with  enhanced  medium  size  particles  is 
accounted  for  by  higher  particle  concentrations  in  the  fine 
and  coarse  particle  range.  This  somewhat  unsatisfactory 
model  will  be  modified  in  the  future  to  allow  for  broader 
size  distributions. 

The  discussion  of  the  details  of  the  lidar  calibration  rou- 
tine shows  that  the  lidar  system  constants  that  have  been 
determined  so  far,  are  valid  only  on  the  very  day  of  their 
determination  if  high  demands  are  placed  on  the  method  of 
remote  aerosol  sensing.  Although  this  is  not  the  final 
stage  of  our  lidar  calibration  procedure,  one  can  at  least 
compute  a sequence  of  aerosol  concentration  and  size  dis- 
tribution profiles  for  a limited  time  interval  of  several 
hours.  Sometimes,  however,  we  found  too  large  variations 
of  the  system  constants  (or  better:  system  factors)  in  the 
course  of  a day  maybe  due  to  thermal  shifts  or  to  almost 
undetectable  damages  in  the  optical  parts  of  the  lidar  sys- 
tem that  are  exposed  to  the  high  energy  of  the  laser  beam. 

In  some  cases  lidar  system  constants  had  to  be  rejected, 
because  vertical  distortions  of  the  backscatter  signals 
led  to  inconsistent  concentration  profiles.  These  vertical 
distortions  may  be  produced  by  poor  adjustment  of  optical 
parts  of  the  system  or  by  considerable  deviations  of  the 
atmospheric  extinction  from  the  values  calculated  after 
Mie  theory.  However,  it  seems  to  be  promising  through  fu- 
ture improvements  at  the  system  to  reduce  the  possible  er- 
rors. 

Some  examples  of  time  sequences  of  aerosol  concentration 
profiles  are  given  in  chapter  IV. 
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2.  Improvements  in  data  processing 

The  evaluation  routine,  which  is  used  to  extract  aerosol 
concentration  and  size  distribution  profiles  from  cali- 
brated lidar  backscattering  signals,  is  not  much  differ- 
ent from  the  calibration  routine  and  was  described  by 
REITER  et  al.  (1975,  1976a).  Some  improvements  of  the 
computer  programs  have  facilitated  the  fitting  of  the 
lidar  data  to  the  impactor  measurements  of  the  aerosol 
size  distribution  at  Wank  peak. 

These  improvements  are; 

- the  automatic  calculation  of  a bimodal  log-normal 
size  distribution  having  the  same  backscattering 
properties  as  a given  size  distribution, 

- an  automatic  fit  of  lidar  backscatter  profiles  at 
694  and  3^7  nm  to  given  backscatter  values,  reduction 
of  the  fitting  parameters  to  the  lidar  system  constants 
and  computation  of  partial  concentration  profiles, 

- an  immediate  graphic  display  of  these  profiles  on 
the  computer  display  terminal  to  check  the  validity 
of  the  lidar  data  and  the  various  parameters  involved 
in  the  evaluation  routine, 

- an  automatic  computation  of  a regression  line  through 
the  logarithmic  lidar  backscatter  profile  in  a select- 
able range  interval  for  comparison  with  Rayleigh  back- 
scatter profiles. 


These  new  facilities  help  to  speed  up  the  throughput  of 
data,  a necessary  condition  for  establishing  time  series 
of  aerosol  concentration  and  size  distribution  profiles. 


IV.  APPLICATIONS  OF  THE  LIDAR  SYSTEM 


1 . Tropospheric  measurements 

The  lidar  system  was  applied  to  the  measurement  of  vertical 
and  time  variations  of  the  aerosol  content  in  an  atmospheric 
column  above  the  lidar  installation.  The  results  of  some  of 
these  measurements  are  presented  in  three  different  ways. 

In  the  case  of  normalized  non-calibrated  backscatter  sig- 
nals, three-dimensional  plots  for  both  lidar  wavelengths 
69A  and  3^7  nm  show  the  qualitiative  variation  of  the  back- 
scattering  medium  with  altitude  and  time.  Examples  are  given 
in  Figs.  IV. 6a,  b - 10a,  b. 

In  some  cases  an  absolute  lidar  calibration  was  carried  out 
and  the  backscatter  signals  were  processed  further  in  order 
to  establish  profiles  of  partial  aerosol  concentrations 
for  fine  aerosol  particles  (about  0.4  ^um  diameter)  and 
coarse  particles  (about  2.0  ^um  diameter)  respectively. 

These  profiles  are  presented  either  in  three-dimensional 
plots  for  both  size  ranges  (Figs.  IV. 4a,  b - 5a,  b)  or, 
for  a better  understanding  of  the  size  distribution,  in 
sequences  of  both  partial  concentration  profiles  in  the 
same  diagram  (Figs.  IV. la,  b - 3a,  b). 


1.1^_  Aerosol  concentration  profiles 

1 ^1 ^1 i Size_distribution_in_|equences_gf _both_particle 

cgncentration_grofileSi 

^§§ic_remarks_on_the_method 

Figs.  IV. la,  b - 3a,  b show  sequences  of  aerosol  concentra- 
tion profiles  in  the  fine  and  coarse  size  range. 

On  8 September  1976  (Fig.  IV. 1)  radiosondes  were  launched 
in  the  morning,  in  the  afternoon  and  in  the  evening  in  or- 
der to  measure  the  changes  of  the  humidity  profile  during 
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the  day.  If  the  relative  humidity  exceeds  70  %,  many  fea- 
tures of  the  lidar  backscatter  can  be  explained  by  particle 
growth  and  changes  in  the  refractive  index  and  the  size 
distribution  due  to  the  water  coating  and  solution  of  the 
hygroscopic  and  soluble  fraction  of  the  aerosol.  In  this 
case  it  is  desirable  to  have  many  humidity  profiles  simul- 
taneously with  the  lidar  firings. 

On  4 March  (Fig.  IV. 2) and  20  May  1977  (Fig.  IV. 3),  the 
aerosol  layers  are  partly  composed  of  Sahara  dust.  Although 
these  cases  are  not  very  appropriate  for  an  absolute  lidar 
calibration,  they  are  good  examples  of  the  applicability 
of  a double  frequency  lidar  to  the  measurement  of  particle 
size  distributions.  Moreover,  the  profiles  on  4 March,  al- 
though rather  "well-behaved”,  show  some  of  the  typical 
features  common  to  all  concentration  profiles  calculated 
with  our  evaluation  routine.  The  features  are  also  inherent 
to  the  size  distribution  model  we  use,  but  to  a high  degree 
resulting  from  a combination  of  the  given  lidar  wavelengths 
and  the  predominant  aerosol  concentrations  and  size  distri- 
butions at  our  measuring  site.  This  will  be  shown  in  the 
following  section. 

Frequently,  the  concentration  profiles  of  fine  and  coarse 
particles  seem  to  be  somewhat  distorted  mirror  images; 
an  increase  of  fine  aerosol  is  accompanied  by  a decreasing 
concentration  of  coarse  particles.  On  4 March  1977.  these 
variations  are  nearly  of  the  same  order  (Fig.  IV. 2a,  b). 

The  two  partial  concentrations  are  not  independent,  they 
are  coupled  through  the  measured  scattering  functions  at 
694  and  347  nm.  Or,  inversely,  the  laser  light  is  back- 
scattered  by  both  fine  and  coarse  particles.  Contrary  to 
the  real  physical  situation,  one  usually  measures  two  in- 
dependent backscattering  functions,  because  the  two  lidar 
wavelengths  use  two  different  paths  through  the  lidar  sys- 
tem (two  independent  calibration  constants)  and  because  li- 
dar firings  in  the  two  wavelengths  are  not  simultaneous. 
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The  partial  concentrations  are  calculated  from  the  follow- 
ing set  of  linear  equations; 


(D  ^ + N2  (D  5 = fj. 
Nl  7=  C 


with 

, N2  : particle  concentration  of  fine  (I)  and 
coarse  (2)  particles, 

; collective  backscattering  function  per  par- 
ticle for  a log-normal  sizo  distribution  with 
mean  diameters  0.4  ^um  (I)  or  2.0  ^um  (2)  at 
wavelength  694  nm  (R)  or  347  nm  (UV) » 

f^’^  : aerosol  backscattering  function  per  unit 

volume  at  wavelength  694  (R)  or  347  nm  (UV)  . 

Collective  backscattering  functions  derived  from  Mie  theory 
are  tabulated  in  REITER  et  al.  (1976b).  For  relative  humidi 
ties  below  50  %,  i.e.  a refractive  index  1.62  - 0.02i, 
these  values  are 

(J)^  = 0.52  10"''^,  <t)^  = 0.67  10"^^, 

(D!j^=  0,26  10"''^,  (J)^=  0,17  10“'^^. 

The  unit  ism  sterad  . 

If  there  is  a range  interval  where  the  backscattering  func- 
tion f?  at  694  nm  remains  nearly  constant  while  the  back- 
" UV 

scattering  function  f^  shows  variations,  the  partial  con- 
centrations and  N2  will  vary  also,  in  opposite  directions 


The  relative  variations  are  of  the  same  order  of  magnitude 
if  the  size  distribution  in  this  range  interval  is  given 
by  the  ratio 


N2/N^  s:  (J)^/  s:  IO'^. 


This  holds  on  March  k,  1977  (Fig.  IV. 2a),  where  the  relative 
variations  of  fine  and  coarse  particle  concentrations  are 
indeed  equal  but  of  opposite  sign. 

The  variation  of  the  particle  concentrations  is  due  to  a 
variation  of  the  backscattering  functions.  As  the  total 
backscattering  function  is  composed  of  aerosol  and  Rayleigh 
scattering,  the  relative  variation  of  the  particle  concen- 
tration may  be  considered  high.  For  instance,  the  relative 
variation  of  the  fine  aerosol  particle  concentration  is  re- 
lated to  the  relative  variation  (or  error)  of  the  backscat- 
tering signal  and  the  system  constant  at  3^7  nm  by 
the  expression 


10®  ( 
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UV 


UV 


.UV 


Q p Tnr 

The  factor  3 . 10  m'^  is  given  by  (J)„  f^  /D,  with  the  Ray- 

IJV  ^ ^ 

leigh  backscattering  function  f^  as  a lower  limit  of  the 
total  backscattering  function  and 

D - ( <0^  0™  ). 


fjjj^  x 8.4  . 10“®  m~^  sterad“\  D s - 1.66  . 10”^®  m^”*  sterad" 

O 

If  the  fine  aerosol  particle  concentration  is  3 . 10  , its 
relative  variation  is  of  the  same  order  as  the  variation  of 
the  backscattering  signal  or  the  error  of  the  calibration 
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constant.  If  this  concentration  is  of  the  order  10  , an 
amplification  by  a factor  of  10  leads  to  large  variations 
in  the  concentration  profile.  This  kind  of  variation  can 
be  found  in  Fig.  IV. la. 

The  considerations  on  error  propagation  do  not  necessarily 
indicate  the  impossibility  of  an  absolute  calibration  of 
the  lidar  system  or  of  the  computation  of  aerosol  concen- 
tration and  size  distribution  profiles  from  calibrated 
lidar  backscattering  signals.  They  simply  show  the  neces- 
sity to  determine  all  parameters  involved  in  the  calibra- 
tion and  evaluation  as  accurately  as  possible,  to  have  a 
most  realistic  scattering  theory  and  model  of  aerosol 
size  distribution.  Furthermore  one  should  be  very  care- 
ful at  interpreting  the  resulting  concentration  profiles. 

1^1 ^2^__Diurnal_va^iations_demonstrated_bY_ three-dimensional 
plots 

Fig.  IV. 4 shows  the  diurnal  variations  as  three-dimensional 
plots  which  were  discussed  before  in  a different  way  in 
chapter  1.1.1.  (Fig.  IV. 2a).  First  a supplementary  remark 
to  Fig.  IV. 2a:  The  lidar  shots  at  11.35  CET  have  been  fit- 
ted to  the  aerosol  particle  data  of  the  station  V ik  ( A ) . 
However,  it  is  satisfying  to  note  that  the  aerosol  concen- 
tration profiles  in  the  coarse  and  fine  range  calculated 
on  this  basis  are  consistent  with  the  aerosol  particle  data 
of  the  stations  Zugspitze  and  Garmisch.  Fig.  IV. 4a,  and  a- 
bove  all  Fig.  IV. 4b  show  in  the  layer  between  2.5  and 
3.0  km  agitations  in  the  aerosol  concentration.  This,  as 
well  as  the  increase  in  the  coarse  range  aerosol  concentra- 
tion with  height  in  the  afternoon  of  4 March  (Fig.  IV. 4b) 
is  indicative  of  heterogeneous  aerosol  conditions  and  of 
an  influx  of  coarser  particles  in  the  height  - a behavior 
being  typical  in  the  case  of  advection  of  Sahara  dust. 
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A similar  behavior  is  particularly  well  demonstrated  by  the 
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fine  and  coarse  particle  profiles  on  22  May  1977  (Fig.  IV. 5) 
when  an  exeptionally  marked  influx  of  Sahara  dust  occurred. 
In  the  course  of  the  day  the  profile  of  the  fine  aerosol 
concentration  (Fig.  IV. 5a)  varies  only  in  the  usual  manner: 
Increase  in  particle  concentration  from  morning  to  after- 
noon in  the  lower  1.5  to  2.0  km,  and  slight  variation  bet- 
ween 2.5  and  3.0  km.  Consequently  the  gradient  of  the  de- 
crease of  the  particle  concentration  with  height  becomes 
increasingly  steeper  in  the  course  of  the  day.  The  profile 
of  the  coarse  aerosol  concentration  shows  a completely  dif- 
ferent pattern  (Fig.  IV. 5b):  Between  1.5  and  2.5  km  we  ob- 
serve a remarkable  rise  with  altitude  the  steepness  of 
which  steadily  decreases  from  morning  to  afternoon,  the 
lower  boundary  of  the  aerosol  layer  with  the  coarse  par- 
ticles continuously  falling  off.  This  behavior  is  quite 
characteristic  of  an  acute  influx  of  coarser  particles  in- 
to the  middle  troposphere.  The  sharpness  of  the  layer  is 
slowly  reduced  by  turbulence  in  the  course  of  several  hours, 
the  lower  boundary  of  the  haze  layer  steadily  decreasing 
until  it  finally  reaches  the  valley  level. 


Just  these  examples  demonstrating  the  influx  of  Sahara  dust 
through  lidar  observations  clearly  show  the  possibilities 
of  applying  a calibrated  lidar  system  in  practice  although 
the  accuracy  presently  achieved  v/ith  regard  to  a quantita- 
tive evaluation  is  not  yet  satisfying  to  a full  extent. 


1 ^2^__piurnal_variation  of  lidar  backscattering  (non- 

three-dimensional 

plots 


Fig.  IV. 6 shows  the  situation  on  20  September  1976  in  ruby 
light  (Fig.  IV. 6a)  and  in  UV-light  (Fig.  IV. 6b). 

We  clearly  observe,  mainly  in  Fig.  6a,  the  increase  cf  back- 
scatter  intensity  in  the  course  of  the  day  in  altitudes  up 
to  about  2.8  km  and  the  rise  of  the  aerosol  from  about  2.0 
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to  2.8  km  and  further  the  formation  of  layers  with  maximum 
backscattering  in  the  range  of  the  upper  boundary  of  the 
exchange  layer  (2.8  km).  This  behavior  is  in  agreement  with 
the  temperature-  and  humidity  lapse  rates  of  the  same  day. 

In  the  afternoon  we  find  the  needle-shaped  peaks  as  a re- 
sult of  precondensation  and  of  the  first,  yet  weakly  de- 
fined cloudiness,  respectively.  The  fluctuations  in  higher 
altitudes  (mainly  in  UV-light)  are  a consequence  of  the  then 
too  unfavorable  signal/noise  ratio  which  quite  recently 
could  essentially  be  improved. Nevertheless,  Fig.  b,  too, 
reveals  (though  somewhat  weaker)  an  increase  of  the  back- 
scatter  intensity  and  rise  of  the  aerosol  in  the  course  of 
the  day. 

Aerosol  structures  can  rather  clearly  be  identified  in 
Figs.  IV. 7 a and  b presenting  the  profiles  of  22  Sep  76. 

This  concerns  on  the  one  hand  the  smoke  swathes  in  the 
forenoon,  but  just  so  the  upper  boundaries  of  haze  remain- 
ing stable  over  longer  time  intervals  with  backscatter  in- 
tensity maxima  in  about  2.1  and  2. 5/2. 6 km  between  14.00 
and  17.00  GET.  These  aerosol  structures  are  only  faintly 
defined  in  UV-light  (Fig.  IV. 7b). 

For  comparison  Fig.  IV. 7c  shows  the  lapse  rates  of  dry  T 
and  wet  T'  bulb  temperature  and  positve  air  conductivity 
for  16.30  GET  on  the  same  day.  According  to  those,  the 
upper  boundary  of  the  exchange  layer  is  in  exact  agreement 
with  the  ruby  laser  backscatter  profile  (Fig.  IV. 7a). 

Fig.  IV. 8 is  also  an  instructive  example  (unfortunately 
without  UV-shot  due  to  malfunction  of  the  instrument). 

It  shows  fairly  well  the  agitation  in  the  aerosol  structure 
up  to  3 km  altitude  as  a function  of  time:  In  the  morning 
a sharp  upper  boundary  of  haze  still  exists  at  2.5  km  which 
is,  however,  soon  being  dissipated.  In  the  valley  floor  we 
have  a wavelike  increase  of  backscatter,  haze  layers  form 
at  various  altitudes  and  the  aerosol  rises  on  the  whole 
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atic  rise  of  the  upper  boundary  of  haze,  and  the  steady 
increase  in  the  backscatter  intensity  at  all  levels  of 
the  haze  layer  is  clearly  to  be  seen.  The  upper  boundary 
remains  sharply  defined  and  at  heights  above  1600  m the 
atmosphere  was  always  free  of  the  haze  advection  (cf.  the 
above  stated  visibility  ranges). 

Fig.  IV. 9c  shows  the  dry  T and  wet  T'  bulb  temperature 
lapse  rates  for  15.30  GET.  The  inversion  lies  between  1150 
and  1300  ra  in  exact  agreement  with  the  lidar  backscatter. 

The  graphical  representation  in  Fig.  9b  was  cut  off  at 
2.2  km  since  the  UV-backscatter  in  the  height  > 2.2  km 
could  not  be  distinguished  from  the  noise  (extremely  low 
aerosol  concentration). 

Anyhow,  it  should  be  pointed  out  that  measurements  of 
aerosol  backscatter  in  ruby  light  are  possible  without 
difficulties  also  under  conditions  of  more  than  100  km 
visibility  as  can  be  seen  from  Fig.  IV. 9a  in  heights 
above  1.6  km. 

We  will  conclude  the  series  of  examples  with  Fig.  IV. 10. 

It  illustrates  fairly  well  the  relatively  rapid  fluctua- 
tions in  the  aerosol  concentration  as  a function  of  time 
and  altitude  without  the  appearence  of  any  stationary  struc- 
tures. The  agitation  in  ruby  light  backscatter  is  conspic- 
uously higher  than  in  UV-light  which  may  lead  to  the  con- 
clusion that  the  inhomogeneities  concern  above  all  the 
coarse  part  of  the  aerosol.  Apart  from  that,  in  the  basic 
structure  good  agreement  is  found  between  Figs.  10  a and 
b.  The  lapse  rates  of  temperature  and  humidity  neither  re- 
veal any  stationary  structures.  This  example  makes  it  clear 
that  a consistent  description  of  the  time-space  aerosol 
structure  can  only  be  obtained  with  the  aid  of  continuous 
lidar  backscatter  profiles  recorded  with  high  time  resolu- 
tion. 
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2.  Observations  of  stratospheric  aerosols  by  photon 
counting  technique 

2^1 Improvements  of  photon  counting  data 

The  off-line  evaluation  of  the  photon  counting  data  remained 
basically  unchanged  since  annual  report  R.  REITER  et  al. 
(1976).  However,  major  improvements  of  the  whole  measuring 
set-up  of  the  lidar  system  are  yielding  much  better  data. 
These  data  are  now  of  such  a quality  that  certain  state- 
ments regarding  aerosol  presence  in  the  lower  stratosphere 
can  be  made. 

Since  October  1976  an  RCA  8852  PMT  can  be  used  instead  of 
the  EMI  9816  PMT  when  utilizing  the  photon  counting  facil- 
ities of  the  lidar  system.  This  new  PMT  is  distinguished 
by  a much  better  time  resolution  resulting  in  a dead-time 
reduction  to  about  one  tenth  of  the  previous  value.  The 
linear  counting  range  now  is  considerably  extended.  With 
regard  to  counting  statistics,  however,  both  PMTs  show 
similar  behavior  which  can  be  expected  within  the  range  of 
total  time  resolution  of  the  single  photon  pulses. 

Another  improvement  is  the  installation  of  a filter  wheel 
in  the  receiver  light  path  holding  several  neutral  filters 
and  providing  unrestrained  light  passage  and  light  attenua- 
tion at  several  levels.  The  filter  wheel  replaced  the  gray 
wedge  which  could  not  be  adjusted  at  reproducible  attenua- 
tion settings.  The  reconstruction  of  the  backscatter  pro- 
file by  computer  calculations  is  very  much  simplified  be- 
cause fixed  attenuation  factors  can  now  be  incorporated  in 
the  program  and  tedious  alignment  of  the  profile  sections 
as  described  in  the  last  report  can  be  ommitted. 

A present  backscatter  profile  is  recorded  by  observing 
600  m intervals  per  channel  of  the  10-channel  counter  and 
7 successive  delay  settings  are  necessary  to  cover  a height 
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interval  from  7 to  40  km.  Per  delay  setting  laser  returns 
are  recorded  until  a 4 96  standard  deviation  is  attained 
but  a maximum  of  50  laser  firings  per  delay  is  not  exceeded. 
Up  to  4 series  per  night  averaged  provide  a mean  profile 
which  consists  of  as  much  as  1000  to  1200  individual  laser 
returns. 

Additional  information  on  the  nature  of  the  scattering 
medium  can  be  gained  by  determining  the  degree  of  depolari- 
zation of  the  return  signal.  An  analyzer  in  the  receiver 
light  path  can  be  adjusted  parallel  or  perpendicular  to 
the  plane  of  polarization  of  the  emitted  laser  pulse,  the 
intensity  ratio  of  the  received  signals  indicating  the 
change  of  polarization. 

2 . 2^__Rayleigh_backscatter  profiles 

The  measured  total  backscatter  (e.g.  molecular  plus  par- 
ticulate backscatter)  has  to  be  compared  with  the  back- 
scatter from  a purely  molecular  atmosphere  (Rayleigh  scat- 
tering). To  minimize  errors  we  try  to  combine  measured 
backscatter  profiles  with  Rayleigh  profiles  derived  from 
actual  atmospheric  density  data  which  are  delivered  by 
our  own  rawinsondes.  These  radiosondes  are  launched  at 
the  site  of  the  institute  at  least  within  12  hours  from 
a lidar  campaign. 

In  case  radiosonde  data  cannot  be  made  available  a 
"standard"  profile  has  to  be  taken  as  a reference.  To 
check  on  the  validity  of  such  a "standard"  profile,  37 
Rayleigh  profiles  obtained  from  radiosonde  data  between 
January  75  and  May  77  have  been  examined.  This  material 
has  been  investigated  regarding  seasonal  trends  and  the 
backscatter  profiles  have  been  classified  as  follows: 
below  the  tropopause  (height  interval  0-8  km)  and  above 
the  tropopause  (height  interval  20  - 30  km)  the  slopes  of 
the  profiles  have  been  determined  and  related  to  the  time 
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scale.  Below  the  tropopause,  which  usually  is  found  bet- 
ween 10  and  13  km,  a seasonal  trend  of  the  inclination  is 
found  with  a maximum  in  June-July  and  a minimum  in  December- 
January.  However,  the  scatter  of  the  data  is  in  the  order 
of  the  seasonal  variations.  Only  little  seasonal  variations 
can  be  found  above  the  tropopause  with  a scatter  by  far  ex- 
ceeding the  variations.  So  any  slope  of  the  profile  can  be 
found  at  any  time. 

Applying  a seasonal  "standard"  Rayleigh  profile  if  actual 
radiosonde  data  are  not  available  therefore  is  only  a 
compromise.  Deviations  of  the  actual  Rayleigh  profile 
from  a "standard"  profile  can  be  in  the  order  of  the 
presently  measured  net  aerosol  backscatter. 

Our  Rayleigh  profiles  have  been  compared  with  the  Rayleigh 
profile  derived  from  density  values  of  the  1976  US  Standard 
Atmosphere.  Below  the  tropopause  there  is  a very  good  agree- 
ment, the"US  Rayleigh  profile"  showing  a slope  corresponding 
to  our  mean  value  of  all  evaluated  profiles.  Above  the  trop- 
opause the  "US  Rayleigh  profile"  shows  a steeper  slope  re- 
sulting in  3.5  % more  molecular  backscatter.  However,  the 
scatter  of  our  profiles  includes  the  slope  of  the  "US  Ray- 
leigh profile"  above  the  tropopause. 


2 . 3^__Measured  backscatter  profiles 

The  quality  of  the  measured  backscatter  data  was  improved 
to  such  an  extent  after  the  introduction  of  the  new  PMT 
that  only  results  obtained  since  October  76  will  be  re- 
ported. Intense  instrumental  maintenance  now  allows  to 
make  use  of  all  clear  sky  periods.  Unfortunately  the  num- 
ber of  nights  suitably  clear  is  rather  limited  and  we  are 
forced  to  conduct  our  observations  during  nights  when  cir- 
rus clouds  tend  to  hamper  our  recordings.  Only  during 
night  time  measurements  can  be  performed  to  avoid  the  high 
background  due  to  diffuse  sky  radiation.  Cirrus  clouds. 
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usually  changing  rapidly  in  density,  can  cause  an  intensity 
shift  of  the  individual  sections  of  the  backscatter  profile 
and  matching  of  these  sections  with  the  general  trend  of 
the  whole  profile  might  then  be  necessary. 

During  the  period  October  76  through  July  77  eleven  back- 
scatter profiles,  in  most  cases  up  to  35  km,  have  been 
recorded,  many  of  them  through  a more  or  less  dense  cirrus 
layer.  These  profiles  are  shown  in  Figs.  IV. 11  - IV. 21. 

The  polygonal  curves  present  the  measured  total  backscatter, 
the  smooth  curves  are  the  calculated  molecular  backscatter. 
The  horizontal  bars  denote  1 standard  deviation  of  the 
statistical  fluctuations  of  the  photon  counts. 

The  matching  of  the  measured  total  backscatter  signal  with 
the  calculated  molecular  return  is  accomplished  at  height 
levels  of  minimum  aerosol  backscatter.  Suitable  levels 
are  above  the  tropopause  around  15  km  and  between  20  and 
25  km.  Deviations  from  the  purely  molecular  backscatter  are 
generally  found  around  the  tropopause  if  cirrus  layers  are 
present,  between  16  and  22  km,  and  above  25  km. 

On  February  14,  1977,  (Fig.  IV. 15)  the  degree  of  depolari- 
zation of  the  return  signal  has  been  determined  between 
5 and  14  km.  The  cirrus  clouds  between  7 and  10  km  could 
be  identified  as  ice  particle  clouds  as  indicated  bv  a de- 
polarization of  60  %.  For  the  molecular  return  a depolari- 
zation of  3.6  % has  been  derived  in  very  good  agreement 
with  published  values. 


2^4^__Stratosgheric_aerosol_ layer 

A good  means  of  describing  particulate  backscatter  is  the 
scattering  ratio,  e.g.  the  quotient  of  measured  total  back- 
scatter to  calculated  molecular  backscatter.  Fig.  IV. 22 
presents  scattering  ratios  derived  from  the  measurements 


between  October  76  and  July  77.  Each  curve  is  the  mean  of 
several  individual  profiles  recorded  during  the  night  of 
observation.  The  tropopause  level  usually  is  found  at  10 
to  12  km  where  in  many  cases  cirrus  clouds  produce  large 
scattering  ratios. 

Above  the  tropopause  a distinct  layer  of  increased  back- 
scatter  intensity  in  the  order  of  10  to  15  % of  the  mole- 
cular return  can  be  observed.  During  the  period  of  obser- 
vation the  lower  boundary  of  this  layer  is  rising  from  a- 
bout  15  to  17  km  and  the  upper  boundary  from  about  21  to 
25  km.  The  maximum  of  backscatter  however  is  reaching  a 
maximum  altitude  at  21  km  in  March  77,  rising  from  18  km 
in  October  76  and  declining  to  19  km  in  July  77.  This 
layer  corresponds  to  that  described  by  Junge  with  scatter- 
ing ratios  reported  for  the  pre-Fuego  situation.  The 
eruption  of  the  volcano  Fuego  in  Guatemala  in  October  197^ 
was  the  last  strong  major  volcanic  injection  into  the 
stratosphere . 

Towards  higher  altitudes  above  25  km  the  scattering  ratio 
is  again  increasing  without  returning  to  unity  at  a de- 
fined altitude.  Although  the  greater  part  of  this  devia- 
tion from  unity  is  covered  by  assuming  an  error  of  10  % 

(to  be  discussed  in  the  following)  our  measurements  seem 
to  indicate  a second  particulate  layer  around  30  km.  Since 
measurements  by  other  groups  at  the  same  latitude  (Garmisch- 
Partenkirchen  47.5°N)  are  not  available,  comparison  and  dis- 
cussion of  our  findings  regarding  a second  layer  are  pre- 
sently not  possible. 


2.5.  Error  discussion 

The  1 0 square  root  error  of  the  photon  counts  amounts  to 
about  -2  96  at  10  to  15  km  and  increases  to  about  -7  96  at 
30  to  35  km.  The  Rayleigh  profile  derived  from  radiosonde 
density  data  will  be  affected  by  errors;  Uncertainties  of 
the  radiosonde  data  will  introduce  an  error,  but  by  matching 


measured  and  calculated  backscatter  profiles  this  error 
will  be  reduced  and  the  remainder  will  not  exceed  1 or 
2 %.  In  some  cases  the  Rayleigh  profile  has  to  be  extra- 
polated to  higher  altitudes  causing  an  error  which  can  be 
regarded  as  of  the  same  order.  For  the  calculation  of  the 
Rayleigh  backscatter  profile  a "no  particulate  extinction 
model"  has  been  applied  which  meets  the  present  strato- 
spheric situation.  The  error  thus  introduced  has  been 
estimated  applying  literature  values  of  the  particulate 
to  extinction  ratio  to  an  assumed  aerosol  layer  between 
18  and  25  km.  An  error  of  .2  % of  the  two-way  transmission 
caused  by  this  layer  has  been  calculated  and  can  be  neg- 
lected. Thus  a maximum  error  of  about  10  % can  be  assumed 
for  our  measurements  between  30  and  35  km. 


V.  CONCLUSIONS 


The  works  under  contract  led  to  the  following  results  or 
conclusions  respectively: 

1 . Computerized,  fully-automatic  multi  frequency  lidar  sys- 
tems make  it  possible,  with  very  high  resolution  regarding 
time  and  distance  (resp.  height)  to  determine  aerosol 
structures  and  to  study  their  variations. 

2.  If  information  is  wanted  relying  not  only  on  relative 
backscatter  data  but  also  comprising  a discrimination 
of  aerosol  concentrations  in  different  particle  size 
ranges,  the  absolute  calibration  of  the  lidar  systems 
is  a necessary  requirement.  This  can  be  done  by  obser- 
vation of  the  pure  molecular  backscatter  as  well  as  by 
in  situ  measurements  of  the  aerosol  size  distribution 
in  the  troposphere  during  lidar  shots.  These  calibra- 
tion procedures  involve  a great  amount  of  effort  and 
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lack  in  their  present  state,  mainly  in  UV-light,  suffi- 
cient accuracy.  Therefore,  further  attention  must  be 
devoted  to  this  subject. 

Nevertheless  some  examples  of  absolute  remote  aerosol 
sensing  show  that  on  principle  the  method  used  is 
practicable  and  delivers  satisfying  results  if  the 
calibration  constants  are  known  with  sufficient 
accuracy. 

3.  Further  it  showed  up  that  even  non-calibrated  lidar 
measurements  in  the  troposphere  yield  valuable  infor- 
mation on  interesting  aerosol  structures  and  their 
variations,  especially  with  regard  to  convection,  ad- 
vection,  and  condensation  processes.  Such  data  can  be 
easily  obtained  with  the  aid  of  fully  developed  auto- 
matic lidar  systems. 

k.  The  method  used  for  single  photon  counting  proved  to 
be  completely  suitable  for  the  performance  of  routine 
monitoring  of  relative  aerosol  backscatter  in  the 
stratosphere  up  to  35  km.  Since  October  1976  we  are 
conducting  at  our  Institute  stratospheric  aerosol 
monitoring.  On  the  average  we  obtain  2-3  profiles 
per  month  during  cloudless  nights. 

The  works  on  the  lidar  systems,  their  calibration  and 
practical  application  will  be  continued  under  contract 
No.  DA-ERO  77-G-064. 
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LEGENDS  OF  FIGURES 

Fig.  II. 1 ; Revised  block  scheme  of  the  stationary  lidar 
system 

Fig.  II. 2 : Block  scheme,  remote  control  electronics 

Fig.  II. 3 : Flow  charts  of  computer  program  for  automatic 

a and  b analog  measurements 

Fig.  II. 4 : Flow  chart  of  the  computer  program  for  auto- 

matic photon  counting  measurements 

Fig.  II. 5 : Laser  with  power  supply  of  stationary  lidar 

system 

Fig.  II. 6 ; Receiver  case  with  radiation  shield  of 
stationary  lidar  system 

Fi.i,  II. 7 : Receiver  case  of  stationary  lidar  system  as 

viewed  from  above,  showing  telescope  mirrors, 
frequency  doubler,  and  45  deg.  mirrors 

Fig.  II. 8 : Data  acquisition  electronics  of  the  stationary 

lidar  system,  including  Biomation  recorder, 
photon  counter,  peak  digitizer,  remote  control 
unit,  and  on-line  computer  with  display  tei'minal, 
hard  copy,  and  floppy  disc  memory 

Fig.  II. 9 : Front  panel  of  remote  control  unit 

Fig.  11.10:  The  mobile  lidar  system 


All  other  figures  are  described  with  full  particulars  in  the 
text . 
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